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ABSTRACT 


This is the second of two papers examining Spitzer Infrared Spectrograph 
(IRS) observations of the ultraluminous X-ray source (ULX) in Holmberg IL Here 
we perform detailed photoionization modeling of the infrared lines. Our analysis 
suggests that the luminosity and morphology of the [0 IV] 25.89 /mi emission line 
is consistent with photoionization by the soft X-ray and far ultraviolet (FUV) 
radiation from the accretion disk of the binary system and inconsistent with 
narrow beaming. We show that the emission nebula is matter-bounded both in 
the line of sight direction and to the east, and probably radiation- bounded to 
the west. A bolometric luminosity in excess of 10 4() erg s™ 1 would be needed to 
produce the measured [O IVj flux. We use modeling and previously published 
studies to conclude that shocks likely contribute very little, if at all, to the high- 
excitation line fluxes observed in the Holmberg II ULX. Additionally, we find 
that the spectral type of the companion star has a surprisingly strong effect on 
the predicted strength of the [O IV] emission. This finding could explain the 
origin of [O IV] in some starburst systems containing black hole binaries. 

Subject headings: black hole physics — galaxies: individual (Holmberg II) 
infrared: ISM - X-rays: binaries 



1. INTRODUCTION 


Despite over two decades of research (Fabbiano 1989), little is known about the intrinsic 
properties of ultraluminous X-ray sources (ULXs). They are characterized by extremely high 
X-ray luminosities, some reaching 10 41 erg s _1 . These luminosities cannot be easily explained 
by normal stellar mass black hole systems, such as those found in our Galaxy. Beamed or 
super-Eddington stellar-mass systems have been proposed as possible models. 
Alternatively, ULXs can harbor intermediate-mass black holes (IMBHs). For a 
complete discussion of various ULX mechanisms, see for example the review pa- 
pers of Miller & Colbert (2004), Soria et al. (2005), and Roberts (2007). Some 
ULX X-ray spectra show distinct features, which have been sometimes inter- 
preted as evidence for the IMBH scenario (Miller et al. 2004). However, the 
super-Eddington mechanism (maybe together with mild beaming) has gained 
ground in more recent years, backed both by theoretical and observational evi- 
dence (Freeland et al. 2006; Done & Kubota 2006; Kuncic et al. 2007; King 2008; 
Soria & Kuncic 2008; Berghea et al. 2008). 

While ULXs a.re easy to detect iu the X-rays they are often faint or undetected in other 
energy bands, making them very difficult to characterize. The ionization nebulae surroun din g 
some ULXs have become critical for understanding the properties of the central source, such 
as the intrinsic luminosity. For instance, Kaaret. et al. (2004) (KWZ hereafter) used the 
optical He II recombination line (A4686) to characterize the intrinsic X-ray luminosities of 
ULXs, as was previously done for a “normal” X-ray binary by Pakull & Angebault (1986). 
Even [Ne V] A3426. a much higher excitation line has been recently detected in a ULX 
(Kaaret & Corbel 2009). However, we note that these and analogous studies are limited to 
the optical regime. 

The very’ luminous ULX in the dwarf galaxy Holmberg II is one such object that has 
been studied extensively in both the X-ray and optical bands (Zezas et al. 1999; Lehmann 
et ai. 2005; Goad et al. 2006; Abolmasov et al. 2007). It is located inside an ionized nebula 
(the “Foot nebula”), and show r s high excitation optical emission lines coincident with the 
X-ray source (Pakull & Mirioni 2002; Abolmasov et al. 2007). The nebula appears to be 
perturbed by a moderate velocity field of ~ 50 km s~ l (Lehmann et al. 2005). The optical 
counterpart of the ULX is a stellar-like bright source located in the “Heel” of the nebula, 
having emission consistent with a B supergiant to a main-sequence O star (Kaaret et al. 
2004). Morever, extended radio emission was detected, approximately coincident with the 
ionized nebula (Miller et al. 2005). The high excitation recombination line He II A4686 
suggests photoionization by the ULX (KWZ), and the radio emission is inconsistent with 
emission from a SNR or an H II region. 



In our first paper Berghea et al. (2009) (Paper I hereafter) we presented Spitzer Infrared 
Spectrograph (IRS) observations of the Foot nebula in Holmberg II. The low resolution 
spectral maps presented in that paper suggest that the high-ionization [O IV] line emission 
is coincident with the ULX position. Indeed we found the mid-IR line ratios for this object 
to be very similar to those seen in standard AGN (Genzel et al. 1998; Sturm et al. 2002; Dale 
et al. 2006). Our preliminary analysis suggested that [O IV], like the He II recombination 
line, can be used to constrain the spectral energy distribution (SED) of the ULX in both 
the UV and X-rays, thereby constraining the intrinsic X-ray luminosity. To this end, we 
presented archival X-ray, ultraviolet (UV) and infrared (IR) data and constructed a detailed 
SED of the ULX from which photoionization models can be used break the degeneracy in 
X-ray models and extract a reliable bolometric luminosity for the source. 

In this second paper we perform detailed photoionization modeling for the Holmberg II 
ULX. We use the spectral energy distribution (SED) obtained in Paper 1 as the ionizing 
source, and compare the predictions with the observed IR Spitzer spectrum. To better 
constrain the physical properties, we investigate how the predictions change when we vary 
the parameters: the accretion disk temperature, gas density, metallicity, and the stellar 
companion. We try both the radiation and matter-bounded geometries to explain the high 
and low-ionization IR. lines. Finally, we explore shock models as a possible alternative to 
ionization. In Sections 5 and 6 we discuss our results and we present our conclusion. 

2. Summary of Observational Results 

The study of ULXs is limited almost exclusively to the X-ray band, where they appear 
very bright and uncontaminated. Studies of these objects at other wavelengths are difficult 
because the emission is generally either undetected or contaminated by circum-stellar or 
circum-nebular emission, making an estimate for the intrinsic bolometric luminosity very 
difficult. This deficiency, especially in the UV, results in degeneracy in X-ray models of the 
accretion process that cannot be broken. Our analysis of the X-ray emission from the ULX 
in Holmberg II is an example of such degeneracy. In Paper 1 three physically distinct models 
were presented to fit the X-ray spectra for the Holmberg II ULX. The resultant parameters 
from the fits are given in Table 1. Figures 1 - 3 show the spectra and the fitted components 
for all three models. A brief summary of the three models is as follows: 

1. Power-law Plus Multicolor Disk (PLMCD) Model: This is a two-component 
model widely used to describe the emission from accreting black holes. As the name 
suggests, the model consists of a multi-color disk component and a power-law compo- 
nent for the hard tail. As discussed in Paper 1 , when extrapolated to softer energies 



(the far UV) this model falls apart since the power law dominates in the UV, Indeed 
the model is unrealistic since it is physically impossible to have the X-ray photons 
Componized to lower energies. This model provides a good fit to the X-ray spectrum 
and is used extensively in other X-ray studies. However we emphasize that it 
is not physically meaningful (Done & Kubota 2006). The PLMCD model is 
shown here as a base from which modified, but more physically meaningful models are 
derived. 


2. The Modified PLMCD Model: The Modified PLMCD model is also a two-component 
model, but one for which the far-UV extrapolation is from the disk component rather 
than the power-law component. This yields a much more physically meaningful inter- 
pretation for the emission. 

3. Broken Power Law (BPL) Model: The spectra were also fit with a broken power- 
law model, which is much more flexible when extrapolated to the UV. Such a model 
would be produced if jet emission dominated the full SED. We note that there is 
currently no direct evidence for jets in ULXs. It is interesting however that 
one of the few sources associated with radio emission, the famous ULX in 
NGC 5408, was shown by Kaaret et ah (2003) to be well fitted by a broken 
PL (see also Soria et al. 2006; Kaaret Corbel 2009). For more details on 
our motivation for these models we refer our reader to Paper 1. 


As can be seen from Table 1, all three models provide reasonably good fits to the X-ray 
data. This degeneracy can only be broken through comprehensive knowledge of the shape 
of the UV continuum and detailed photoionization modeling of high-excitation lines. In the 
former case, the three X-ray models can be extrapolated to the UV and then compared 
with observational data. The best-fit extrapolated model to the UV spectrum provides 
some indication of the appropriate physical accretion scenario for the observational data. 
However, contamination to the UV emission by circum- stellar or circum-nebular continuum 
precludes a pure fit. Thus, in addition to this photometric data, high-excitation spectral 
lines, which require high-energy UV photons to ionize, can be used in combination with 
photoionization modeling to check the robustness of the fit. In Paper 1 we also presented 
multi-wavelength photometric data from archival Spitzer . XMM -Newton and GALEX ob- 
servations of the Holmberg II ULX, and also Hubble Space Telescope and Very Large Array 
photometric data from the literature. These data were used to construct a detailed spectral 
energy distribution (SED) of the ULX (Fig. 4). All three X-ray models enumerated above are 
shown in the plot as dotted, thick solid, and dashed black lines respectively. The thick solid 
black line in Figure 4 represents the "base SED model” which is the input spectrum for the 
preliminary photoionization modeling explored below. As discussed in Paper 1 and shown 



in Figure 4. the Modified PLMCD model is chosen as the X-ray model for the base model, 
since it is the only extrapolated model that fits the optical and UV data well. Tables 1 and 2 
provide the X-ray spectral fitting results and the IR line measurements from the first paper 
that will be the basis of the modeling in the current publication. 


3- CLOUDY Modeling of the IR lines 

To perform detailed photoionization modeling with CLOUDY 1 , a number of input pa- 
rameters need to be defined; Specifically, the geometry of the cloud needs to be described, 
the density and metallicity of the gas needs to be specified, and an input SED needs to be 
established. The following outlines the specific input parameters used for the preliminary 
CLOUDY run. These parameters are later varied to check the final results. 

L Gas Cloud Geometry; The input geometry must describe the shape and thickness 
of the gas cloud. Shape includes, for example, a spherical, cylindrical, or plane-parallel 
geometry. The cloud thickness defines a radiation or matter bounded geometry. In 
their photoionization modeling, KWZ assumed a spherical geometry based on resolved 
optical lines, and assumed the nebula is radiation-bounded. However, they noted that 
the lower-ionization line [O I] is only detected west of the ULX, thereby suggesting a 
matter-bounded nebula in other directions. Though we also detected low-ionization 
lines in the IR spectrum (Table 2 in this paper and Figs. I and 2 in Paper 1), we 
begin with the spherical, radiation-bounded geometry in an attempt to reproduce 
the high-excitation |0 IV] emission. We assume the inner radius of the cloud to be 
0.1 pc, thus maximizing the ionizing flux. To constrain the size of the ionized nebula 
for the radiation bounded geometry, CLOUDY was run down to a temperature of 
4000 K. At this temperature, the cloud thickness is 140 pc, and the column density is 
4.2 xlO 21 cm™ 2 , more than ten times the X-ray estimated value. 

2. Density: Density can be specified as varying, uniform or clumpy. In the preliminary 
CLOUDY model, we assume a constant gas density of 10 cm™ 3 and a filling 
factor of 1 as per KWZ. This density (and filling factor) was estimated by KWZ 
using the optical line surface brightness profiles predicted by CLOUDY. 

3. Metallicity: Metallicity is expressed in terms of solar metallicities. We chose an input 
metallicity of 0.1 Z@, close to that of KWZ and Pilyugin et al. (2004). 


1 Photoionization modeling was performed using version 07.02.01 of CLOUDY, described by Ferland et 

al. 1998. 



4. SED: The input SED must specify the companion star and the X-ray model used to 
fit the UV to X-ray continuum. The input SED used in our preliminary analysis is 
described in Section 2 and results from the multi- wavelength photometric and spec- 
troscopic observations presented in Paper 1. This so called ’’base SED model” utilizes 
the Modified PLMCD model, extrapolated to the UV. For the input SED we used a 
B2Ib supergiant spectrum to represent the stellar companion. The B2lb supergiant 
was chosen because it better fits our UV data than the 05V (both are limits to the 
stellar spectral range suggested by KWZ). The bolometric luminosity of this model is 
1.34 xlO 40 erg s” 1 . 


3*1* Results from the Base Model 

The results of the first CLOUDY run, with the base models as the ionizing source, are 
given in Figure 5 (top- left panel), as the ratio of the predictions to the observed values. The 
CLOUDY model predicted values for the mid- IE lines are close to the observed values for the 
high ionization tines. A simple change to the base model such as a slightly younger companion 
star would be enough to fit the [O IV] line observed, as the analysis below indicates. We 
also show predictions for the He II flux measured by KWZ (2.7 x 10 3b erg s _L ). 

The low ionization lines are the most discrepant. This disparity between the model and 
the observations is likely caused by a matter-bounded geometry in at least some directions 
of the nebula, as already suggested by KWZ. However, many of the model parameters are 
not well known and we explore how changes to these can effect the line-flux predictions in 
the sections to follow. 

We note that sensitivity and aperture effects could explain the discrepancy between the 
low ionization line fluxes predicted by the radiation-bounded base SED model and the ob- 
served values. Using the CLOUDY output we calculated surface brightness profiles expected 
for the mid-lR lines. CLOUDY predicts a spatial extent of 20 to 25 pc for the He II which 
is consistent with the findings of KWZ (15—26 pc). The [O IV] emission from the ULX 
should be extended to a radius of about —20 pc according to the simulation. However, due 
to the poor resolution of our spectral maps (one LH pixel = 67.5 pc), it is unclear whether 
the model accurately predicts the spatial extent of the [O IV]. In addition, the CLOUDY 
models suggest that while sensitivity may have some effect on the observed low- ionization 
line fluxes, it is not solely responsible for the high predicted-to-observed line ratios for the 
radiation bounded base SED model. 



3.2. Varying the X-ray Source Spectrum to the BPL and PLMCD models 


We next tried changing the X-ray part of the source spectrum. We first replace the 
Modified PLMCD model with the PLMCD model described in Section 2, leaving the other 
CLOUDY parameters the same. The model is extrapolated to the edge of the [O IV] line 
where we introduce a break, so that it does not dominate the optical flux from the B2 star. 
Results from this CLOUDY simulation are shown in upper-right panel of Figure 5. The 
values predicted by this model are clearly not consistent with the measured values. This 
model predicts a strong [Ne V] line, which is not detected in our Spitzer observations, and 
a He II luminosity 27 times larger than measured by KWZ. The [0 IV] flux predicted from 
this model is also an order of magnitude larger than the observed flux. We therefore rule 
out this model as a good estimate for the X-ray source spectrum. We note that the PLMCD 
model has been used by many researchers to estimate disk temperatures and even black hole 
masses. Our findings imply that it is likely not an appropriate model for ULX high-energy 
spectra. 

We then tried replacing the Modified PLMCD model with the BPL model described in 
Section 2. Here we found results similar to the PLMCD model, with [Ne V] predicted to be 
detected, and all the observed lines well overpredicted. We therefore conclude that neither 
the PLMCD model nor the BPL model with a radiation bounded geometry provide a good 
fit to the observed IR spectra. This was expected because both these models predict a very- 
strong UV ionizing flux. 


3.3. Varying Other Parameters 

We investigate here in more detail how the strong radiation from the accretion disk 
photoionizes the surrounding ISM and produces the observed lines. Starting again with the 
base SED model in the radiation bounded geometry, we vary the input parameters to see 
how they affect the emitted spectrum. We first examine the dependence of the simulated line 
luminosities on the X-ray spectrum by changing the MCD disk temperature (kT iTt ), then the 
gas density, metallicity and the inner radius of the cloud. Finally, we explore how the SED 
of the stellar companion affects the line luminosities. The models presented are based on the 
multi-wavelength data of the Holmberg II ULX and the known properties of the surrounding 
ISM, but this analysis can be applied to ULXs in general. 

Varying the Disk Temperature: The disk temperature is first altered by shifting 
the X-ray base SED model to lower photon energies, such that the disk luminosity remains 
constant across the band. As can be seen from Figure 6a, when the disk component peaks 
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at lower energies, the mid-IR emission lines become stronger. For cooler disks (0.1 keV) the 
[Ne Vj line should be detected in our observations. In contrast, the lower ionization lines 
are not dependent on the disk temperature over the range usually found in ULXs (0. 1—0.5 
keV) . We conclude from this that varying the disk temperature will not yield better results. 

Varying the Gas Density of the Cloud: The results obtained by varying the density 
of the cloud are plotted in Figure 6b. The line luminosities are not particularly sensitive to gas 
density (see also KWZ), though the [O IV] luminosity increases slightly as the gas approaches 
the critical density. KWZ found that the He II line is clearly extended to 15 -“25 pc and 
consistent with a density of ~10 cm” 3 . We therefore conclude that varying the density of 
the gas will not do much to improve the model. 

Varying the Metallicity: The dependence of the line luminosities on metallicity is 
shown in Figure 7a. The [O IV] and [Ne V] line dependence is linear. Interestingly, the base 
model with a higher input metallicity reproduces the measured [0 IV] luminosity very well. 
In fact a metallicity of 0.19 Z 0 (oxygen abundance of 9.3xl0~ 5 in absolute value) reproduces 
the [O IV] to within the calibration error of the instrument, but without increasing the He II. 
which is already over-predicted in the current version of the model. We conclude from this 
that the input metallicity may be lower than it should be (see Pilyugin et at. 2004). However, 
we also note that a higher metallicity will predict higher fluxes for the low ionization lines, 
resulting in an even greater discrepancy for these lines. We also explored the effects of 
gas phase depletion on metallicity and found that depletion had very little effect 
on either the metallicity or the resultant luminosity of the [O IV]. 

Varying the Inner Radius of the Ionizing Cloud: Figure 7b shows that the ionizing 
power of the incident flux is diluted if the spherical cloud is too far away from the source. 
If the inner radius is larger than about 10—20 pc, the luminosities of the [O IV] and [Ne V) 
lines drop dramatically. Interestingly, the He II line, though similar to [O IV] in ionization 
potential, is nearly insensitive to the inner cloud radius. This is likely because the [O IV] and 
[Ne Vj are in a higher stage of ionization than He II and as a result are much more affected 
by dilution of the ionizing photons (in other words a drop in the ionization parameter) from 
increasing the inner radius. We conclude from this that increasing the inner radius will only 
worsen our original results and will not affect the strength of the low ionization lines. Thus 
we conclude that the inner radius is likely close to the source as our original input parameter 
supposes. 

Varying the Stellar Companion (the 05V star): To test the effect of the stellar 
companion on the IR lines we ran the base model with an 05V star replacing the B2Ib 
supergiant. The temperature and luminosity of the 05V star, 42000 K and 3.2 xlO 39 erg s -i , 
were chosen to match the model used by KWZ. The results (the lower- right panel of Fig. 5) 
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show that the 05V star has a very strong impact on the high ionization lines. The luminosity 
of the [O IV] line increases by a factor of 3, and the [Ne V] line by factor of 2. If the companion 
is in fact an 05V star, our measured Optical Monitor (OM) and GALBX fluxes from Paper 1 
are under-luminous and are likely affected by reddening. Inclusion of an 05V star also has 
a very different effect on the [O IV] versus the He II line. Indeed, the He II luminosity is 
predicted to increase by only 11%. The strong impact of the 05V star can be explained if the 
O star produces copious 0 ++ ions (i.p. 35 eV), which are in turn ionized further by the ULX. 
The O star thus acts as a catalyst to the higher-energy' photons from the ULX. A comparison 
between the oxygen ions structure for these models is presented in Figure 8a. The number of 
photons in the 24.6— 54.4 eV interval is three times higher for the 05V star (7.27xl0 48 s~‘ 
compared to 2.42 xlO 48 s -1 ) than the B2Ib supergiant. Figure 8b show's that the local cloud 
temperature for the 05V star is also about two times higher in the region where the 0 3+ 
ions are produced (at offsets of 10—20 pc). Indeed, the [O III] A5007 luminosity estimated 
from CLOUDY increases by a factor of six when the 05V star is used instead of the B2Ib 
supergiant. 

The companion star therefore has an indirect but strong impact on the [0 IV] line 
emission. A similar effect can be seen on a smaller scale for the [Ne V] line, which has an 
even higher ionization potential (97 eV). Thus, the [O IV] line does not act simply as a high 
energy photon counter as is the case with the He II line, but can also provide information 
about abundances of lower energy photons. This effect can be used to explain the detection 
of [O IV] in star-forming regions. In this scenario, [O IV] can be produced by relatively 
faint sources of soft X-rays and UV photons hidden in the star forming region, including 
such sources as ULXs, X-ray binaries, or SNR.. This is a significant result, since the origin 
of [O IV] in star forming galaxies has been the subject of debate in the literature (Lutz et 
al. 1998; Schaerer & Stasinska 1999). 

Because our CLOUDY modeling overpredicts most of the high ionization 
line fluxes when an OV-5 star is used, we conclude that an OV5 star is likely not 
the companion star for the Holmberg II ULX. On the other hand, later-type 
stars, such as an 08 star (or later), might better predict the [O IV] flux in the 
base model. In support of this, if the ULX is powered by a stellar-mass BH, the 
companion required to reproduce the observed X-ray luminosities must be very 
massive (Podsiadlowski et al. 2003). Thus a stellar mass black hole could explain 
both the [O IV] luminosity and the very high X-ray luminosity. The fact that 
the Holmberg II ULX resides just outside a UV~bright region of star formation 
(see Paper 1), also suggests that the environment is well suited for formation 
of such massive stars (see also NGC4559 X7 Soria 2007). However, without the 
black hole mass constrained, an IMBH with a massive companion is still also a 



plausible option. 


3.4* Varying the Geometry 

As mentioned at the beginning of this section, KWZ found evidence from low- ionization 
optical lines that the ionized nebula is matter-bounded, at least in the east. This is also 
suggested by the IR line maps presented in Paper 1 and by the over- prediction of the low- 
ionization lines by the radiation bounded base SED model above. Paper 1 also shows that 
the low ionization line fluxes are likely contaminated by circum- stellar emission* 
We also noticed that the absorption suggested by the X-ray spectrum is too low compared 
to that implied by a radiation-bounded geometry. More precisely, for our base SED model, 
the simulation requires a column density of 4,2 xlO 21 cm" 2 , which is more than ten times 
larger than the measured values from spectral fits (see also the beginning of this section). 
This suggests that the nebula is also matter- bounded in the line-of-sight direction, 

To simulate a mat ter- bounded geometry, we first ran the base SED model with the 
geometry described by a spherical cloud of varying thickness, so that the total hydrogen 
column ranges between 10 20 and 10 22 cm" 2 . These results are plotted in Figure 9a. The 
column obtained from the X-ray fit is shown as a vertical line and labeled “Matter Bounded” . 
For this value, the predicted fluxes for both [O IV] and He II lines are much lower than the 
observed values. Therefore, this simple matter bounded geometry" using the base SED model 
is not sufficient to reproduce the observed lines. As Figure 9a shows, the minimum column 
density required for the base model to predict the measured values for the high excitation 
lines is ~5 xlO 20 , or a cloud thickness of 16 pc (the dashed vertical line). This 16 pc is 
within the 15-25 pc radii observed in the He II emission line maps from KWZ* Therefore the 
matter bounded case could work, if the observed X-ray column density were slightly higher 
than the X-ray predicted value of 3.26 xlO 20 cm -2 . 

We next replaced the Modified PLMCD model for the BPL model as the X-ray fit for the 
SED. The broken power-law model is better suited than the Modified PLMCD model to test 
the matter bounded geometry. By varying this absorption parameter in the model within 
the errors allowed by the fit, we obtained a family of broken-power models, all consistent 
with the X-ray data, but predicting very different slopes in the ionizing part of the SED 
(the photon index varies between 1*3 and 2.2). A subset of the broken power-law models 
are shown in Figure 4. We re-ran our CLOUDY simulation for these models with the 
corresponding column densities from the X-ray fits. The results are plotted in Figure 9b* 
The best fit broken power-law X-ray model (BPL) is shown as a vertical line in this figure. 
The predicted line luminosities for this model are given in the lower-left panel of Figure 5. 



As was the case with the original, radiation bounded base SED model, the predictions 
are very good for the high-ionization lines. The cloud radius is 17.7 pc, which is consistent, 
with the size of the nebula found by KWZ in the east and south. However, contrary to the 
original radiation-bounded base model, the ULX makes no contribution to the low-ionization 
lines. This is true even for the largest column density, 9.2 xIO 20 cm -2 (a 29.7 pc cloud radius 
for a gas density of 10 cm” 3 ). However, the spectral maps in Figure 3 of Paper 1 show that 
the ULX is likely to contribute at least to [Ne III] and probably [Si III. 

Thus we conclude that the matter- bounded geometry alone does not fit all of the mid-IR 
spectral observations well. In fact, these results suggest that the geometry is complicated, 
being asymmetric and including both radiation and matter bounded geometries. KWZ 
found similar results. They find that to the south and east, He II and H/3 emission suggests 
a matter- bounded geometry: the low excitation lines are absent here but detected iri the 
west at larger offsets. Our CUBISM maps from Paper 1 suggest something similar, though 
the spatial resolution here is limited compared with the HST observation. The fact that 
the low-ionization line luminosities predicted by the matter-bounded models are much lower 
than the ones predicted by the radiation-bounded geometry suggests that a combination of 
the two geometries would yield predicted values that are very close to the observations. 


4, Shock Models 

The radiation-bounded / matter-bounded geometric scenario presented in Section 3 to 
explain the mid-IR line fluxes in the Holmberg II ULX presumes photoionization as the source 
of the line emission. However, shocks could also produce such lines. Indeed, a combination 
of X-ray photoionization and shocks was recently proposed to explain the morphology of 
optical lines observed in the vicinity of LMC X-T (Cooke et al. 2007). In this case the shocks 
are likely driven by a jet. Abolmasov et al. (2007) performed optical spectroscopy on eight 
ULXs with nebular counterparts, including the Holmberg II ULX, and found evidence for 
shocks in all cases based on the [S II] AA671 7.6731/Ha ratio. However, using high-excitation 
lines such as He II and [O III], these authors also found that Holmberg II ULX and three 
other ULXs required photoionization to sufficiently explain the optical ratios. 

Lehmann et al. (2005) found evidence of radial velocity variations of ± 50 km s” 1 from a 
region similar of size to the He III nebula. It is therefore possible that part of the ionization 
is caused by shocks, Lutz et al (1998) suggested that velocity shocks of 1 00—200 km s'™ 1 can 
produce significant [O IV] emission. However, when we also take into account the optical 
data, and in particular the detection of the He II recombination line, our data provide strong 
evidence against shocks as the source of the [O IV] emission. This is primarily because the 
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He II is not easily produced by shocks. Instead, the simultaneous detection of these lines 
argues for the presence of significant radiation over the He II edge (see also Schaerer & 
Stasiriska 1999), and therefore suggests photoionization as the most likely line production 
mechanism. 

Nevertheless, we checked to see if shock models could explain the mid-IR line luminosi- 
ties and ratios observed in Holmberg II. We used Mappings III grids from Allen et al. (2008) 
to construct diagnostic diagrams. The predicted [() IV] emission from shocks is shown in 
Figure 10a. We show both pure shocks and shocks with ionized precursors at a gas density 
of 1 and metallicity 0.1 Z®. For these models, the parameters varied are the shock 

velocity and the “magnetic parameter” B/n 1/2 , where B is the magnetic field and n is the 
gas density (Dopita & Sutherland 1996). 

The CUBISM measured fluxes from the brightest and faintest pixels in the [() IVj map 
(yellow and blue pixels in Fig. 3 of Paper 1, respectively) are shown as horizontal lines in 
Figure 10a. We note that the maximum intensity of the line is likely higher than these 
levels. Indeed the (O IV] emission is not spatially resolved and the plotted levels are actually 
averages over the size of the LH pixel (67.5 pc), rather than the maximum from each of the 
averaged slits. In spite of this, the levels plotted in Figure 10a requires significant shock 
velocities and such velocities (> 200 km s _1 ) were not detected by Lehmann et al. (2005). 

For optical lines we chose the He II and the [O III] AA4363,5007 lines, as per Evans 
et al (1999). The ratio [O III] A4363 / [O III] A5007 is sensitive to the gas temperature. 
Using the data from Lehmann et al. (2005) we calculated the ratio ([() III] A4959 + [O III] 
A5007) / [O III] A4363 = 50, and estimate a temperature of 17800 K (Osterbrock & Ferland 
2006). This is consistent with the temperature of the gas in the region of the cloud where 
0 ++ ions are produced (see Fig. 8b). For the He lI/H/3 ratio we chose an average (0.2) 
between the measurements of KWZ, Pakull & Mirioni (2002), and Lehmann et al. (2005). 
The shock grids for the optical lines are plotted in Figure 10b. For comparison, we over 
plotted the published data for two other ULXs with detected nebulae, Holmberg IX ULX 
and MF16 in NGC 6946 (Abolmasov et al. 2007). While the Holmberg IX ULX data is 
consistent with shocks with modest velocities (< 100 km s -1 ), Holmberg II ULX and MF16 
seem to require very high velocities. Even for the shock plus precursor models, velocities in 
excess of 300 km s -1 are required to reproduce the observed line ratios and Lehmann et al. 
(2005) find velocities of only 50 km s -1 . 

In summary, some contribution from shocks to the high excitation lines detected from 
Holmberg II ULX cannot be ruled out completely. However, if there is some contribution 
from shocks to the high ionization emission line luminosities, it is likely very small, since 
both He II and jO IV] lines require shocks with velocities in excess of 200 km s~ l , and such 
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velocities were not seen by Lehmann et al. (2005). 


5. Discussion of the Modeling Results 

We find that the Spitzes observations of the Holmberg II ULX and especially the detec- 
tion of the (O IV] tine is consistent with photoionization by radiation from the ULX. Of all 

the lines examined, only the [O IV] was shown in Paper 1 to be correlated with 
the ULX. However, we also found in Paper 1 that some contribution by the ULX 
to the [Si II], [S III], and [Ne III] emission line flux is also likely. The modeling 
results presented here are based on the observations presented in Paper 1 and 
are summarized as follows. 

1. The Base Model (Modified PLMCD): The CLOUDY simulations show that the 
Spitzes spectrum is consistent with photoionization by the accretion powered emission 
from the ULX. We find that limited sensitivity might explain the relatively low' fluxes 
observed for some low-ionization lines. The base SED model predicts a slightly lou r er 
[O IV] luminosity than the measured values for the radiation bounded case (Fig. 5), 
and even low'er for the matter-bounded case (Figure 9a). Previous X-ray results have 
shown that the accretion disk temperature is likely lower than in our base model (KWZ 
obtained ~0.2 keV) . Such a disk temperature wall increase the predicted line luminosity 
to match our measurements (Fig. 6a). A slightly higher metallicity (0.2 Z©) can yield 
better fits to both the [O IV] and He II data for this base model. This is because 
the He II line is already over-predicted and raising the metallicity will not increase 
this discrepancy (Fig. 7a). In support of this, a metallicity of 0.17 was measured by 
Pilyugin et al. (2004). We note that while this model reproduces the spectra very well, 
it requires absorbing columns in the X-rays that are inconsistent with the observations. 

2. The PLMCD Model: We find that the widely used PLMCD model, w r hen extrap- 
olated to the UV, is not consistent with the observations. The high excitation lines 
predicted using this model are a.t least ten times larger than observed. In this case 
lowering the metallicity to match the [O IV] will not work, because He II will still be 
much over-predicted. This model has been used extensively to fit the X-ray data from 
ULXs (e.g. “cool disks” Miller et al. 2004). However, the X-ray community has 
more recently begun to fit these spectra with more physically motivated 
models Freeland et al. (2006); Done & Kubota (2006); Kuncic et al. (2007); 
Soria (2007); Soria U Kuncic (2008, e.g.), which appear to be consistent 
with high quality ULX X-ray data (Berghea et al. 2008). That the PLMCD 
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model is inconsistent with the Spitzer observations of the Holmberg II ULX 
supports these previous X-ray studies. 

3. Matter vs. Radiation Bounded Geometry: The intensity and morphology of 
the 1R lines provides spatial information about the emission line nebula around the 
Holmberg II ULX, The geometry consistent with our infrared data (and previously 
published optical data) is approximately spherical but not symmetrical. The CUBISM 
maps show that to the east and south the nebula shows only high excitation litres 
that axe likely matter-bounded, while in the west, the detection of lower excitation 
lines might suggest a radiation-bounded cloud. However, observational effects, such as 
aperture and sensitivity effects, preclude any definitive conclusion on this model. The 
discrepancies found in Section 3.1 between the CLOUDY model and the measured 
values for the low- excitation lines suggest that the matter-bounded geometry is at 
play at least in part (probably to east and south as found from the spectral maps). 
Moreover, the absorption measurements from the X-ray data require the cloud to be 
at least matter-bounded toward our line of sight. Therefore, we find that both matter 
bounded and radiation bounded geometries are likely at work to describe the observed 
spectra in both the Mid-IR and X-rays. 

4. Shock Models: The shocks models for both IR and optical emission lines show that 
shocks with velocity >200 km s'" 1 are required to generate both the He II and [O IV] 
lines observed in Holmberg IL The detected radial velocity variations around the ULX 
are of the order of ±50 km s^ 1 (Lehmann et al. 2005). Contribution from shocks to 
these lines is therefore likely to be small. 


5. True X-ray Luminosity: The ULX luminosity estimated using photoionization mod- 
eling is independent of the estimate based on the X-ray flux, which can be affected 
by both absorption and beaming (King et al. 2001). The nebular emission line lu- 
minosities, on the other hand, are time-averaged true luminosities and can provide 
information about the geometry of emission and the surrounding ISM. 

It is now widely accepted that ULXs are variable on timescales from days to years. 
These timescales are much shorter than the ionization equilibrium timescales in the 
nebula, Chiang & Rappaport (1996) investigated time-dependent photoionization of H 
and He for supersoft X-ray sources. The He II ionization equilibrium at the 90% level is 
reached in ~0.7 r R after source turn-on, where r R — 3000 yr is the recombination time 
for He ++ from KWZ. For periodic sources, it was found that the He II line luminosity 
decreases significantly with increasing source period. For example, for a period of 
10 r R and a duty cycle of 10%, the estimated luminosity is 4.7% of that expected from 
a steady source with the same peak luminosity. That is a 53% deviation from a linear 



scaling law (10% of the steady state). This implies that if the Holxnberg II ULX is 
variable, the estimates of KWZ based on the He II line are actually lower limits of the 
time- averaged lumi nosity. 

For 0 + '\ the recombination time is shorter, r# ~ 100 yr (Osterbrock & Ferland 
2006), for the same ionization conditions: electron density 10 cm" 3 and temperature 
20000 K. We expect therefore that the [0 IV] line reaches equilibrium in a shorter 
time (< 100 yr), but still long enough to provide a time-averaged estimate of the true 
luminosity of the ULX over ~10 - 100 yr. We note that light travel time (53 yr for a 
0 +a nebula of 16 pc), is an important factor in this case, because it is comparable to 
the recombination time. 

Finally, the morphology of the emission lines, from both the previously published 
optical data and our IR data suggest emission that is nearly isotropic, being inconsistent 
with narrow beaming. Moreover, photoionization modeling with our base model is 
consistent with the detected high ionization lines, the true bolometric luminosity of 
this model is >10 40 erg s"U This is well above the Eddington luminosity for a stellar- 
mass BH. In addition, some observational results indicate that the BH in Holmberg II 
is likely much larger than a stellar mass black hole (see the conclusion below). 

6. Impact on IR Studies of Star bursts: We find that the predicted [O IV] line 
luminosity is significantly affected by the type of the companion star. Faint [O IV] 

. has been detected in many star bursts and star-forming regions, but its origin is still 
under debate. Our analysis provides a new mechanism for this emission, suggesting 
that the [O IV] line could be produced by relatively faint X-ray sources including X-ray 
binaries, SNR or ULXs in starburst galaxies. 


6 + Conclusion 

The bolometric luminosity of the base model is 1.33 x 1G 40 erg s”U For accretion within 
the Eddington limit, the mass of the central BH is at least 85 The estimate based on 
the inner disk temperature (0.38 keV) in the MCD model gives a much larger mass, 994 M 0 . 
Finally, scaling from AGN estimates, the [O IV] emission predicts a mass as large as 10 4 M 0 . 

All of the estimates presented here presume the black hole is accreting within 
the Eddington limit. However, while our analysis excludes strong beaming in 
Holmberg II ULX, a sub-Eddington IMBH or a super-Eddington stellar-mass BH 
are both plausible ionization mechanisms that stem from our analysis. Indeed 
we show in Section 3.3 that a massive companion (e.g. B2Ih star) and a stellar 
mass black hole can reproduce the [OIV] luminosities see in the Spitzer spectra, 



as can an IMBH and a similar companion. The CLOUDY modeling results alone 
cannot discriminate between these two sources or the two black hole estimates, since we have 
no way of constraining the accretion rate. 
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Table 1. X-ray model fits 
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MOD flux 
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Modified PLMCD 

3.26*0,56 

0.38*0-02 

2.54*0.02 

1.02/135, 04)74/713 

1.11 

0.49 

PLMCD 

10.4±0,70 

0.27*0.02 
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•2.28 

0.12 

BPL 

5-4T.s-0.7B 

1.73*0. 10 

2. 57 ±0.02 

1.0316/649 

1.4.1 

U-60 


Note. (H: X-ray model, as described in Section. 2. (2); Intrinsic hydrogen column density. The Galactic 
column density from KWZ (3.42 x It? 20 cm 2 ) was added separately. (3): inner disk temperature for the MOD 
component. For the broken power-taw model (BPL), this is the first (low energies) photon index parameter. 
The break is at 10 keV. (4): Photon index for the power-law component. For the broken power-law model this 
is the second (high energies) photon index parameter. (5): Reduced >f 2 values for the fit and the number of 
degrees of freedom. For the Modified PL XI CD model, we quote two valuta for each separately fitted component 
(see text for details). (6); U nub sorbed (intrinsic) luminosities between 01 and if) koV. (7): XI CD component 
e.ui absorbed fiux as fraction of the total fiux for the PLMCD model- For the other two models, this simply the 
ratio of the flux between 0.1 and 1.0 kcV to the (lux in the whole range (0.1 -1.0 keV) 


Table 2. Measured infrared lines 
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Note. We show line fluxes* foe the standard (slightly smaller than the SH map) mid the 
jen&s) aperture (4 LH pixels) as defined in Section 2 of Paper 1. Fluxes are in 113 W cm 
luminosities (L) in 10 3 ^ erg s"' 1 . If the measurements errors are smaller than the absolute 
calibration accuracy of 25%, the latter were used. For uon.de t notions we show 3c upper limits. 
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PI-MCD Model 



Photon energy (keV) 


Fig. 1. PLMCD model fit to the XMM-Newton spectrum. This is the typical model used 
for ULXs: a multi-co 1 -- — 1 * 


Modified PLMCD Model 



Photon energy (keV) 


Fig. 2. — The Modified PLMCD model fit. Here the MOD component and the power-law 
are fitted separately below and above 1 keV, respectively. 
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BPL Model 
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Fig. 3. — The broken power- law model fit. The break is at 1 keV, the same energy where 
the two components axe separated in the Modified PLMCD model in Fig, 2, 
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Fig. 4. The input SED (Base model) constructed in Paper 1 and used as ionizing source 
in CLOUDY. The radio data is from Miller et aL (2005), the V-band magnitude for the 
optical counterpart of the ULX is quoted from KWZ. Our measured fluxes in Paper 1 using 
IRAC, MIPS, OM and GALEX are upper hounds, and therefore our model is below these 
data points. The two stellar spectra correspond to the spectral types range consistent with 
the colors and magnitudes in KWZ. The three models fit to the X-ray data are shown, 
extrapolated to the [O IV] edge at 54.93 eV (shown as a vertical line). The BPL model 
was used to obtain a family of broken power-law models, all compatible with the X-ray data 
within the errors. We show the lower and upper limits of this family as the lower and upper 
dashed lines, corresponding to F -- L3 and F -- 2,2 in the low energy slope, respectively. 
The middle dashed line represents the best fit BPL model shown in Table 1, The SED is 
constructed using the prefered X-ray model. Modified PLMCD. 
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Ratios Predictions / Observations for Various Input SEDs 



Fig. 5. We show predictions for the IR lines from the CLOUDY modeling for various 
ionizing SEDs. We also included the He II A4686 line measured by KWZ (2.7 x IQ 36 erg s^ 1 ). 
They are presented as model-to-data flux ratios. For non- detected lines, the plotted values 
are lower limits and are marked with an arrow pointed upward. The predictions for the 
[Ne III] and [Si II] lines are < 0.01 in the BPL model (lower-left) and are not shown. The 
perfect match line (ratio — 1) is shown on ail plots. The errors are small, similar in size 
to the data points, therefore are not shown. The measured line fluxes are presented in the 
Table 2. 
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MCD kT a (kc V ) |j>g Gas Density (e*n ) 

Fig. 6. The dependence of the simulated lines from CLOUDY on: a) MCD inner disk 
temperature, b) Gas density of the cloud. The vertical lines mark the base model values: 
kT m -■ 0.38 keV, density = 10 cur 3 . The line colors are the same but only labeled in the 
left panel 
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Fig\ 7. The dependence of the simulated lines from CLOUDY on: a) Cloud rnetallicity. b) 
The inner radius of the cloud. The vertical line marks the base model value of Z — 0.1 Z 0 . 
The line colors are the same but only labeled in the left panel. 
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Fig. 8.- Comparison between the base model, with the B2lb star (continuous line) and the 
model with the 05V star (dashed line), a) The oxygen ionic structure in the cloud (0 2+ and 
0 3 b). b) Local temperature in the cloud. 
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Fig. 9. This figure illustrates the difference between a matter bounded and radiation 

bounded nebula. The line colors are only labeled in the right panel, a) The dependence of 
the simulated lines from CLOUDY on the cloud hydrogen column density for the base model. 
The equivalent depth of the spherical cloud is also marked on the top X-axis. Both vertical 
continuous lines mark the column density for the base model, but one value was obtained by 
fitting the X-ray data with the base model (matter bounded nebula), while the other is the 
column in the CLOUDY 7 simulation where the local temperature reached 4000 K (radiation 
bounded). We also marked with a vertical dashed line the limit where the predictions for the 
high excitation lines are still compatible with observations, b) The simulated lines using the 
family of broken power-law models fitted to the X-ray data. The best fit broken PL model 
(BPL) is marked by the vertical line. The bottom X-axis shows the low-energy photon index 
for the models. On the top X-axis we mark the corresponding log column densities (lower 
values) and cloud thickness (upper values). All these models are matter-bounded. 
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Fig. 10. Mappings III shocks and ionization models. The shock grids are plotted for 

magnetic parameters: 0.5, 1, 2, 3.23, 5 and 10 // G cm 3 / 2 and for shock velocities 125, 250, 
300, 400, 500, 700 and 1000 km s -1 . a) [O IV] line intensity shocks predictions. The 
shocks plus precursor models vary little with the magnetic parameter, and we only show the 
model with a nominal equi partition value (3.23 fiG cm 3 / 2 , the dashed line). The horizontal 
lines correspond to the brightest (continuous) and respectively faintest (dotted) pixels in 
the [0 IV] line map in Figure 3 of Paper L b) Optical lines diagnostic with shocks and 
photoionization models. The pure shock models are in red, shocks plus precursor in blue 
and the photoionization grid in black. For the photoionization grid we varied the density 
(1, 10 and 100 cm -3 ), and the inner disk temperature (0.1, 0.2, 0.4, 0.8 and 1.0 keV), and 
assumed a metallicity of 0.1 Z 0 . For Holmberg IX ULX and MF16 we show with arrows 
the approximate locations for Solar metallicity. While Holmberg IX ULX is consistent with 
shocks with modest velocities (< 100 km s -1 ), Holmberg II ULX and MF16 seem to require 
very high velocities. They are however consistent with photoionization. 





